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Abstract 


We apply heavy-quark effective theory (HQET) to separate long- and short- 
distance effects of heavy quarks in lattice gauge theory. In this paper we focus 
on flavor-changing currents that mediate transitions from one heavy flavor 
to another. We stress differences in the formalism for heavy-light currents, 
which are discussed in a companion paper, showing how HQET provides a 
systematic matching procedure. We obtain one-loop results for the matching 
factors of lattice currents, needed for heavy-quark phenomenology, such as the 
calculation of zero-recoil form factors for the semileptonic decays B —> D^lv. 
Results for the Brodsky-Lepage-Mackenzie scale q* are also given. 
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I. INTRODUCTION 


This paper applies heavy-quark effective theory (HQET) to study the renormalization 
in lattice gauge theory of currents containing heavy quarks. It is a sequel to our papers on 
power corrections [l]] and on radiative corrections to heavy-light currents ||. Here we treat 
the case where one heavy quark flavor decays to another. To make contact with heavy- 
quark phenomenology, we denote the two flavors by the labels b and c. In the description 
of lattice gauge theory with HQET, discretization effects of the heavy quarks are absorbed 
into the short-distance coefficients of the effective Lagrangian and effective currents. The 
key difference between this work and Ref. J2[] is that now, with two heavy flavors, HQET is 
used to describe both heavy quarks. Thus, the short-distance coefficients are functions of 
m;,a and m c a, where mb and m c are the heavy quark masses, and a is the lattice spacing. 

Heavy-heavy currents are needed to calculate the zero-recoil form factors of the semi- 
leptonic decays B —> Dlv and B —> D*lv, as well as the change in the form factors away from 
zero recoil. These decays are of great phenomenological interest, because with reliable lattice 
calculations of B —> transitions one can make a model-independent determination of 

the element V c b of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. For further background 
on the impact of lattice matrix elements on CKM phenomenology see, for example, Refs. ||. 

The application of HQET to lattice gauge theory began with the consideration of 
power corrections |lj], building on the demonstration [|J that Wilson’s formulation of lattice 
fermions [^] have a well-defined heavy quark limit. In particular, the Isgur-Wise heavy quark 
symmetries 0 emerge whenever the heavy-quark mass Aqcd, even if mn ~ a -1 . (The 

label h denotes a generic heavy flavor.) Because Wilson fermions have the the same sym¬ 
metries and fields as continuum heavy quarks, on-shell lattice correlation functions can be 
described with HQET. This point of view is sometimes called the “non-relativistic interpre¬ 
tation” of lattice QCD ||. 

In describing either lattice gauge theory or continuum QCD with HQET, physics at short 
distances is lumped into short-distance coefficients. Several short distances—two inverse 
heavy quark masses and (on the lattice) the lattice spacing—arise in the decay of one heavy 
flavor to another. The coefficients depend not only on and m c a, but also on other 
“irrelevant” parameters of the lattice action and currents. By adjusting these parameters 
of the action and currents, lattice gauge theory can be tuned, term by term, to the heavy- 
quark expansion of continuum QCD |4].|T]||]. In this adjustment, details of HQET, such 
as its renormalization scheme, drop out. Matching with HQET is, thus, an intermediate 
conceptual step that explains how to match lattice gauge theory to QCD when m^a ^ 1. It 
provides an attractive extension of more familiar matching procedures, such as those based 
on the Symanzik effective field theory [pj-jTTI 


which usually assume m^a <C 1. 


In addition to developing the formalism, which holds to all orders in perturbation theory, 
we compute the one-loop terms of the matching coefficients of the leading dimension vector 
and axial vector currents. We present compact expressions of the integrands of the momen¬ 
tum integration, for the Fermilab action and currents [Q. We also present numerical results 
for the Sheikholeslami-Wohlert (SW) (11|] and Wilson || actions, which are the special cases 
most often used in the (non-relativistic interpretation of the) Fermilab heavy-quark method. 
We include the so-called rotation terms [|T^,[4| in the currents, which are needed for tree- 
level matching at dimension four. It is also possible to obtain most of the normalization 
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non-perturbatively [[L3|-[l5|], and we give separately the residual short-distance correction. 


15 


16 


Our perturbative results have been used in the calculation of the zero-recoil form factor 
for B —> D*lu [L| and, indirectly, in the calculation for form factors for B —> nlu 
One-loop results without the rotation were given in a preliminary report of this work 
and used in the calculation of the zero-recoil form factors for B —> Dlv p~3 . 

This paper is organized along the lines Ref. Section |J discusses how to use HQET 
to separate long- and short-distance physics with (continuum) effective field theories. In 
particular, we obtain a definition of the matching factors of the heavy-heavy vector and 
axial vector currents. Section E introduces a specific definition of lattice currents suited 
to improvement in the HQET matching procedure; these currents generalize those used 
recently for B —> D ( *' } matrix elements |L|,[l4j. In Sec. [TV] we present one-loop results 
for the matching factors, including the scale q* in the Brodsky-Lepage-Mackenzie (BLM) 
scale-setting prescription JT^|T^|. Some concluding remarks are made in Sec. 0. Three 
appendices contain details of the one-loop calculation, including an outline of a method to 
obtain compact expressions, and explicit results for the one-loop Feynman integrands for 
the vertex functions with the Fcrmilab action. 


Instead of printing tables of the numerical results in Sec. IV. we are making a suite of 

nrsn 


programs freely available 
treated in our companion paper 


This suite includes programs for the heavy-light currents 


II. MATCHING WITH HQET 


In our companion paper [2] we reviewed how the standard Symanzik description of cutoff 
effects breaks down when m^a ^ 1. It is worth re-emphasizing that it is the description 
that breaks down—particularly the Taylor series of short-distance coefficients in powers of 
m^a. Lattice gauge theory remains well-defined for all m^a, but one needs other tools to 
understand how to relate lattice observables to continuum QCD. The obvious alternatives 
are the effective field theories HQET pH -24] and NRQCD |25|,26], which exploit the simpler 
dynamics of systems with one or more heavy quarks. The simpler dynamics also emerge in 
lattice gauge theory with massive fermions, so the effective theories also can be re-applied 
to understand lattice observables Ell¬ 
in this section we show, in the case of heavy-heavy currents, how to use HQET to match 
lattice gauge theory to continuum QCD. We first recall the HQET description of heavy-heavy 
currents in continuum QCD, paralleling the discussion in Ref. |2|. We then explain what 
changes are needed to describe lattice gauge theory with heavy fermions. We focus on HQET 
because of the phenomenological importance of B —> D W transitions; for quarkonium the 
logic could be repeated with NRQCD velocity-counting rules |25| , |26|] . LInlike the standard 
Symanzik program [0 0 , the HQET approach works even in the region where m^a is no 


longer small. Like the usual HQET, however, it requires (p is a typical momentum) 


m h >p, Aqcd, (2.1) 

but once this condition holds (and pa ■ Cl), our formalism provides a systematic description 
of lattice observables for all m^a, h = 6, c. 
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The conventions for HQET are the same as those given Sec. Ill of Ref. [|f]. Each HQET 
quark field carries a velocity label. For the time being we use two different velocities, v 
and v', for the two flavors. The velocities can be chosen somewhat arbitrarily, but HQET 
is a good description of QCD if each is close to the velocity of the hadron containing the 
heavy quark.[| The fourth Euclidean component V\ = iv°, so in the rest frame v = (i, 0 ), 
and similarly for v'. The metric is taken to be diag(±l, 1,1,1), with the upper (lower) sign 
for Euclidean (Minkowski) spacetime. In either case, v 2 = v' 2 = —1. We denote the two 
HQET fields as c v > and b v . When the flavor of quark is unimportant, we write formulas with 
the more generic symbol h v . The HQET field h v satisfies the constraint |(1 — if)h v = h v , 
or 

Tfb v = ib v , Crftf = ic v ,. (2.2) 

Physically this constraint means that only quarks, but not anti-quarks, are described. The 
tensor r)% — 8% + v p ‘v u projects onto components orthogonal to v. For a vector p, the 
component orthogonal to v is pj_ = rjtfp 1 ' = + v ,J 'v ■ p. In the rest frame, these are the 

spatial components. Similarly, r/'() = 6% + and p^_, = rj'^p 1 ’. 

HQET describes the dynamics of heavy-light bound states with an effective Lagrangian 
built from c v > and b v . For each flavor one writes 

£qcd = ~ q(P + m)q = ^hqet, (2.3) 

where the symbol = means “has the same on-shell matrix elements as”. The HQET La¬ 
grangian consists of a series of sets of terms 

£hqet = £ (0) + £ (1) + £ (2) + • • •, (2.4) 

where consists of all operators of dimension s + 4 built out of h v and h v (and gluons 
and light quarks). The ultraviolet regulator and renormalization scheme of the two sides of 
Eq. ( |2.3|) need not be the same. For the aims of this paper it is enough to consider the first 
two terms, and £ (1 ) ; but the generalization to higher dimension should be clear. 

The leading, dimension-four term is 

P 0 ' 1 = h v (iv ■ D — mh)h v , (2.5) 

where rn^ is the mass of flavor h. is a good starting point for the heavy-quark expan¬ 
sion, which treats the higher-dimension operators as small. The mass term in £ ( °) is usually 
omitted, because it obscures the heavy-quark flavor symmetry (though only slightly [|l|). By 
heavy-quark symmetry, it has an effect neither on bound-state wave functions nor, conse¬ 
quently, on matrix elements. It does affect the mass spectrum, but only additively. 

When the mass term is included, higher-dimension operators are constructed with = 
D M — iMv 11 and V = D — iMv’. Here M selects the mass of flavor h: Mh v = m^hy. 


1 In NRQCD, the relative velocity between the heavy quark and heavy anti-quark of quarkonium 
should not be confused with the velocities v and v 1 introduced here. Note that it is also possible 
to formulate NRQCD with a total velocity label like v 


25], 
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To describe on-shell matrix elements one may omit operators that vanish by the equations 
of motion, —iv ■ T>b v = 0 and c v 'iv' • V — 0, In practice, higher dimension operators are 
constructed with V ^ = Dj_ (or, for velocity v ', Dj_,) and [D>\ T> ,J ] = [D^,D U ] = F^ u . The 
dimension-five terms in the Lagrangian are 

£d) =C 2 0 2 +C B 0 B , (2.6) 

where C 2 and Cb are short-distance coefficients, and 

0 2 = h v Dlh v , (2.7) 

Ob = h v s aP B a0 h v , (2.8) 


with s a /3 = —ia a p/2 and B a ^ = T^rpF ,w . In operators with two HQET fields (for flavor¬ 
changing currents, say) T>' and T> appear together. 

I 11 Eq. ( |2.5|) one should think of the quark mass rrih as a short-distance coefficient. By 
reparametrization invariance |29| and Lorentz invariance of (continuum) QCD, the same 


mass appears in the short-distance coefficient of the kinetic energy C 2 O 2 , namely, 


C 2 = 


2 rrih 


(2,9) 


If the operators of HQET are renormalized with a minimal subtraction in dimensional reg¬ 
ularization, then rrih is the (perturbative) pole mass. The chromomagnetic operator Ob 
depends on the renormalization point of HQET, and that dependence is canceled by 




2 m h ' 


( 2 . 10 ) 


The description of flavor-changing currents proceeds along similar lines. Let 


= cip l b 


( 2 . 11 ) 


be the flavor-changing vector current, where c and b without the velocity subscripts are QCD 


fields obeying the Dirac equation. In HQET is described by 

14 

= CV||U M c<A, + C v± c v 4^b v + C Vv ,v'lc v fb v - ^ B Vi Q^ i -|-, (2.12) 

2=1 

where c v > and b v are HQET helds, which satisfy Eq. (|2.4|) and whose dynamics are given 
by £hqet- The fourteen dimension-four operators are 

Qvi = ~v»c v ,p L b v , (2.13) 

Qv 2 = c V ’i r y±P±b v , (2.14) 

Qv 3 = c v 'iDj_b v , (2.15) 

Qv 4 = +v^c v >p ±l b v , (2.16) 

Qv 5 = c v/ p ± ,H±fb v , (2.17) 
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Qve = c v >iD^,b v , (2.18) 

0.V7 = -v'^c v 'p ± b v , (2.19) 

Qvs = ~ v±cvP ± ,b v , ( 2 . 20 ) 

Qv 9 = —v^c v iiv' ■ D ± b v , (2.21) 

Qv io = c v 'i^j_iv' ■ D±b v , ( 2 . 22 ) 

Qyn = -v^c v 'iv ■ D±fb v , (2.23) 

Qv 12 = D/w • -DjR7±A, (2.24) 

Qvi 3 = —v^Cu'iv’ ■ D ± b v , (2.25) 

Qyi 4 = -v'£c v 'iv ■ D±ib v . (2.26) 


Further dimension-four operators are omitted, because they are linear combinations of those 
listed and others that vanish by the equations of motion. In developing the heavy-quark 
expansion the appropriate equations of motion for the heavy quarks are (—iv ■T>)b v = 0 and 
c v r(iv' ■ V) = 0, derived from the respective £ (0) . 

The QCD axial vector current is 

A» = 6 , (2.27) 

where c and b are again Dirac fields. A has an HQET description analogous to Eq. (|2.12|) , 

14 

A 11 = C A± c v 'i'Y±'Y5b v ~ CA^Cy’^by - C Av ,v'lc v ^b v -• ( 2 -28) 

Z— 1 

By convention, each operator Q!' m is obtained from Qy i by replacing c v > with — c^/ 75 . 

The coefficients Cj and Bji and the operators Qji play an analogous role to the coeffi¬ 
cients Cj and Bji and the operators Qji introduced in Ref. [2J, but they are not the same. 
The latter are defined in an effective theory using a Dirac field q (without a velocity label) 
to describe the light (er) quark, whereas the barred symbols are defined in a theory using an 
HQET held c„' (with a velocity label) to describe the lighter (heavy) quark. The bars are 
used here to emphasize the difference. 

There are many fewer operators when v' = v. One set of operators vanishes because 
v'j_ —* v± = 0, namely, Qj 7 , Qj S , Q!j Vi - Qju, another set vanishes because v' -D±, v-D±> —> 
v ■ Dj_ = 0, namely, Qj-, j > 9. A last set vanishes because P + (v)TP + (v) = 0, where T 
stands for the full Dirac structure sandwiched by c v and b v . In the end, one is left with only 
v^Cyby, Qy 2 , Qy 3 , Qy 5 , and Qy 6 for the vector current, and c v i^j 5 b v , Q^, and Q% for 
the axial vector current. 

The coefficients of the HQET operators depend on the heavy-quark masses and m c 
(to balance dimensions), as well as m c /rrib and /i/mj, where /i is the renormalization scale. 
At the two-loop level and beyond, dependence on light quark masses rn q also appears. With 
two velocities, the coefficients also depend on w = —v' ■ v. Although they are not needed 
below, it is instructive to give the coefficients of the dimension-three terms. At the tree level 
= C'j] = 1 and C j 0 ', = 0. Through one loop order, for v' — v, 
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(2.29) 


Cvj, = 1 + C F ^^3f(m c /m b ), 

C A± = 1 + C F ^l [3 f(m c /m b ) - 2], (2.30) 

where 

f(z) = ^jlnz-2. (2.31) 


The important properties of f(z) are /( 1) = 0, f(l/z) = f(z). At this order, /i dependence 
appears only in the renormalized coupling g 2 (/i). The “non-/” part of C A± given here 
assumes that the renormalized axial current satisfies the axial Ward identity [[30]. At the 
tree level, the coefficients of the dimension-four currents are 


dM _ dM _ ^ _ l dM _ l dM /'O 

t>V\ — &V2 — — +-R41 — +W42 (Z.6Z) 

d [°] _ d [ 0 ] _ ^ _ _ d [°1 _ _ S [°1 (2 QQ ') 

n V4 — n V5 — 2 m — U A4 — \z.66) 

B [ j] = 0, otherwise, (2.34) 


but all become non-trivial when radiative corrections are included. 

HQET provides a systematic way to separate the short distances 1 jm c and 1 /rrib from the 
long distance Aq CD , as long as the condition ( |2.1| ) holds for both flavors. HQET can also be 
applied to lattice gauge theory with the same structure and logic, treating a as another short 
distance. The strategy is viable as long as condition ( | 2 . 1 | ) holds and, of course, pa - Cl. For 
lattice NRQCD applied to heavy-light systems, HQET is just the corresponding Symanzik 
effective held theory. HQET may also be applied to Wilson fermions, because they have the 
correct particle content and heavy-quark symmetries as QCD [|I|]. Thus, for a lattice gauge 
theory with either NRQCD or Wilson quarks 


Aat — ^HQET- 


(2.35) 


The effective Lagrangian £hqet has the same long-distance operators as in Eq. (|2.3| ). The 
lattice modifies the short-distance behavior, so lattice artifacts of the heavy quarks should 
be lumped into the HQET coefficients [|IJ]. For example, in the dimension-four HQET La¬ 
grangian £(°) 

, one must replace rrih with the lattice rest mass m^. 

Starting with £ (3) one must introduce operators to describe lattice violations of rota¬ 
tional invariance. An example is the dimension-seven operator O 4 = h^(D/) 4 /i^. Such 

operators do not, of course, appear in the HQET description of continuum QCD, but they 
do not upset the general framework. They can still be defined as insertions in a continuum- 
regulated theory (analogously to operators like I] At ?(D M ) 4 g in a Symanzik effective held 
theory). Because the symmetry-breaking arises at short-distances, it is more useful to focus 
on the operators’ coefficients, noting that C 1 ^ 7 ^ 0, instead of C 4 = 0 as for continuum QCD. 

Returning to £ (0 ) and A 1 / the coefficient of the kinetic energy becomes 


C lat 
2 


1 

2 m 2 h 


(2.36) 
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and rri 2 h is called the kinetic mass. If the HQET operators are defined by minimal subtraction 
in dimensional regularization, then both and m 2 /, can be computed from the pole in the 
perturbative quark propagator |33| , and they are infrared finite and gauge independent 
The lattice breaks Lorentz (or Euclidean) invariance, so reparametrization invariance no 
longer requires m- 2 h to equal mu,. 

For the flavor-changing currents, one forms bilinears of lattice fermions fields with the 
right quantum numbers. Any such lattice currents for b —> c can be described by 



14 


V& = C^v^b,, + C%c v 'i'y1b v + Cy*v±c v 'b v - £ + • • •, 

i —1 

Kx = c%_csi-rZ'rsK - - £ sSai + ■ • • 


i =1 


(2.37) 

(2.38) 


An explicit construction of V{,^ t and A Xai through dimension 4 is given for Wilson fermions 
in Sec. [TTJ, and f° r lattice NRQCD already in Ref. ||32|| . On the right-hand side the HQET 
operators Qj % are the same as in Eqs. ( |2.12|) and ( |2.28| ), but the short-distance coefficients 
Cj at ^ t and B l jl are modified. In particular, they depend on the lattice spacing a through m^a 
or m c a, and also on adjustable improvement parameters in V^ x and A XaX . At and beyond the 
dimension six, there are HQET operators to describe violations of rotational invariance in 
the lattice currents, but, as above, they do not spoil the overall framework. 

Since Eqs. (|2.13|) - (|2.26|) give a complete set of dimension-four HQET currents, the coef¬ 
ficients (7j at J± and B ] y contain short-distance effects from both heavy quarks. For HQET 
operators of dimension s + 3, the corresponding coefficient must contain s powers of the 
short distances m^ 1 , m~ l or a. Since they are functions of all ratios of short distances, it 
is a matter of choice which dimensionful quantity is factored out. The key point is that, 
even when applied to lattice gauge theory, the heavy-quark expansion is still a successive 
approximation: higher terms in the expansion are suppressed by powers of Aqcd, with the 
dimensions balanced by the short distances. 

Lattice artifacts from gluons and light quarks can be described by the Symanzik formal¬ 
ism M3 . At some level the light quarks will influence the HQET coefficients of the heavy 
quarks, and the heavy quarks will influence the Symanzik coefficients of the gluons and light 
quarks. These are, however, details that do not obstruct the central idea of using effective 
field theories to separate long- and short-distance physics. 

The abstract ideas can be made more concrete by comparing the HQET descriptions of 
continuum and lattice matrix elements. The continuum matrix element of v ■ V is 


(D\vV\B) 


= -Cv,{D 

flf) - V B Vi (D^\v e«|s<°>) 


*eS|| 

— C2cCv }] j 

’ d 4 x{D^\T0 2c (x)c v/ b v \B^y 

— CbcCv I, J 

f d A x(D^\TO Bc {x)c v ,b v \B^y 

— C 2 bCy |! J 

" d 4 x(D^\Tc v/ b v 02 b (x)\B^y 

— C-BbCy || J 

f d A x{D^\Tc v ,b v O Bb {x)\BW Y + 0(A 2 /m 2 
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(2.39) 
















where the set S n = {1,4, 9, 11}, and B (D) is any fr-flavored (charmed) hadronic state. (The 
*-ed T product is defined in Ref. JjJ; this detail is unimportant here.) On the left-hand side 
the states are QCD eigenstates; on the right-hand side (with velocity subscripts) they are 
the corresponding eigenstates of £ (n \ the leading HQET Lagrangian. The extra subscripts 
on the coefficients and operators from C < ' ]} denote flavor. The corresponding lattice matrix 
element is similar, but slightly different: 


Viatl B) = —Cy*{D\ 

*S5|| 

/■7 lat/Mat i 
U 2c °V]| I 

f d^x(D^ \T0 2c (x)c v ,b v \BW)* 

/Mat/Mat i 
( M3c°Vj| I 

( d A x(D^\TO Bc {x)c v ,b v \BWy 

/^lat/^lat i 
U 2b / 

f d 4 x(D^\Tc v/ b v 0 2b (x)\B^y 

/Mat/Mat i 

| J 

f d 4 x(D^\Tc vl b v O Bb (x)\B^y 

+ K a . F C ^ 

Jd 4 x(D^ \Tc v ib v qiaFq(x)\B^y + 0(A 2 a 2 b(ma )) 


(2.40) 


Now the states on the left-hand side are lattice eigenstates. But on the right-hand side 
the differences compared to Eq. (|2.39 ) are all in the coefficients: except for the last matrix 
element, they are, term by term, the same. The last T-product arise from the Symanzik 
local effective Lagrangian for light quarks, cf. Ref. ||. 

Similar formulas hold for matrix elements of V F and and for the axial vector 

current. If one multiplies the lattice matrix elements with 


z h 


c, 


/'-flat ’ 
°Jn 


c 


J ± 


/-flat ’ 


(2.41) 

(2.42) 


and subtracts the lattice from the continuum equations, a simple picture of cutoff effects 
emerges: lattice artifacts of the heavy quark are isolated into the mismatch of the short- 
distance coefficients, namely, 


lat 


sc t = a 

SB Jt = ZjiB% 


c 

i ■ 


B 


Ji- 


(2.43) 

(2.44) 


In the lattice term of 5Bji, matching factors Zj t appear to restore a canonical normalization 
to the lattice currents. One has Zji = Zj for i e S\\, and Zji = Z J± for i G S± = 
{2,3,5,6,7,8,10,12,13,14}. 

The matching factors Zj^ and Zj ± play the following role, sketched in Fig. [Ij. The 
denominator converts a lattice-regulated scheme to a renormalized HQET scheme, and the 
numerator converts the latter to a renormalized (continuum) QCD scheme. As long as the 
same HQET scheme is used, HQET drops out of the calculation of Zj n and Zj . Moreover, 
changes in continuum renormalization conventions modify only the numerator, and changes 
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lattice 


C/C ut 

QCD 



FIG. 1. Diagram illustrating how the matching factors C' lat , C, and Z 
gauge theory and QCD to HQET, and to each other. 


C/C lat match lattice 


in the lattice action or currents modify only the denominator. Thus, the matching factors Zj n 
and Zj ± play the same role as the factors Zj^ and Zj ± in HQET matching for heavy-light 
currents 0. 

When HQET matching is applied to Wilson quarks, it is possible that the lighter heavy 
quark satisfies m c a -C 1. Then one could equally well apply the heavy-“light” formalism 
of Ref. 0, describing the charmed quark a la Symanzik with a Dirac field. This regime 
is interesting because in practice one often has m c a < | and A/2 m c a < ~, so both the 
Symanzik and the HQET descriptions are reasonably accurate. 

By comparing the two matching procedures, one can derive relations between the two 
sets of short-distance coefficients. 2 To proceed, one must introduce coefficients to relate the 
heavy-“light” theory with the Dirac field c to the heavy-“heavy” theory with the HQET 
field c,/: 


cb v = Gy, c v >b v + Y B Vi v ■ Q Vi H-, (2.45) 

ieS'n 

ci'y±b v = C Vl c v >i r i±b v + C Vv ,vfc v >b v - Y ^vi VvQvi 4 -> (2.46) 

ieSx 

14 

Qvi = Y, CvijQvj + • • •, * = 1,6, (2.47) 

3 = 1 


and similarly for the axial vector current. The relation between the three sets of short- 
distance coefficients —C and B. C and B, and C and B —is sketched in Fig. [^. By substi¬ 
tuting Eqs. ( |2.45|) -( f2~47l ) into Eq. (2.40) of Ref. @ and comparing to the above Eq. ( j2.37|) , 
one hnds (when m c a <C 1) 


/^lat /^lat ry 

Wj L “ ( - y J ± ( ^J± 

/olat sy\a,t sy 

L J V , ~ c Jx C A' 


(2.48) 

(2.49) 

(2.50) 


2 One must take the same lattice currents, 
nearly the same. 


The improved currents of Ref. ||2[ and Sec. Ill are 
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lattice 



lat 


QCD 

cTb 



FIG. 2. Diagrams illustrating how the matching factors C, C, and C match the underlying 
theory to effective theories where the lighter quark is described either with a Dirac field c or an 
HQET field c v > (note subscript in HQET case); for (a) lattice gauge theory, (b) continuum QCD. 


B'j‘< = E B'gCjij + C'fBjj, j e S,, 


2=1 

6 


Bfi = E B'jICjh + C'f’Bjj, j 6 Sj 


(2.51) 

(2.52) 


2=1 


The same relations hold for the coefficients describing continuum QCD, i.e., dropping the 
superscript “lat”. The same coefficients C and B appear, because no lattice fields appear in 
Eqs. ( j2.45| )- (PT~nj ). In Fig. g, the bottom sides of the triangles are oblivious to the underlying 
theory, be it lattice gauge theory or continuum QCD. Thus, one can eliminate the Cj from 
Eqs. (|2.48|) - (|2.50|) in favor of C/C, yielding relations between the matching factors, when 


m c a -C 1, 


Z J\\ - z J\v 

z j± = 

ry /-ylat 

c J v /• 


(2.53) 

(2.54) 

(2.55) 


In particular, Eqs. (|2.53|) and (|2.54|) show that the heavy-heavy matching factors and Zj ± 
are a natural extension, into the regime m c a ^ 1, of the heavy-light matching factors ■Z'y 
and Zj ± . Consequently, in the following sections we shall drop the distinction and omit the 
bar for the heavy-heavy matching factors. 

One also finds relations among the coefficients of dimension-four currents: 


6 


- Bjj = E (Zj,B% ~ B M )C Mj , j 6 S lh 

2=1 

(2.56) 

ZjJjI - Bjj = jz(Zj L B'}l - j 6 S ± . 

2=1 

(2.57) 


The terms on the right-hand sides are proportional to mismatches in the heavy-light HQET. 
Each ZjB 1 ^ — Bji is suppressed by a power of a, multiplied by a function of m b a. Similar 
constraints can be derived for higher-dimension coefficients. Following this reasoning for the 
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whole tower of higher dimension operators leads to the conclusion that, when m c a -C 1, the 
mismatch of the heavy-“heavy” currents reduces to that of the heavy-“light” case. In this 
way, lattice gauge theory reproduces the entire l/m c expansion, apart from deviations that 
are suppressed by powers of the other short distances, a and l/m*,. 

Deep in the a —> 0 limit, where m c a < m^a <C 1 one can relate the short-distance 
coefficients of the HQET formalism to those in the Symanzik formalism.[] For simplicity 
let us take v = v'. For the [axial] vector current, one inserts Eq. ( [2.12D [Eq. (|2.28| )] into 
Eq. (2.5) [Eq. (2.6)] of Ref. 0, neglects terms of order m 2 a 2 , and compares with Eq. (|2.37|) 
[Eq. (|2.38|)]. One also must match the tensor bilinear to HQET at the dimension-three level, 


cia^b = C T r]^c v ia a \ 


+ * • * 


(2.58) 


where one introduces another short-distance coefficient Ct■ At the tree level, C y = 1. The 
pseudoscalar bilinear is described by HQET operators of dimension four and higher, because 
P + (u) 75 P + (u) = 0. After carrying out these steps, one Ends that 

£vj, = Zvi 

Z V \\ By\ = Byi + aZ v K v C T , i = 2, 6, 

Zv n Bv\ = B V3 — aZyRyCx, i — 3, 5, 

from matching the vector current (with Zy and K v as dehned in Ref. |2|), and 

= Z A , (2.62) 

Za ± B 1 21 — B A i + 0(a~), i = 1,4, (2.63) 

from matching the axial vector current (with Za as dehned in Ref. [2j). Consequently, one 
sees 

lim Zj^ = Bji, (2.64) 

a —>0 


(2.59) 

(2.60) 
(2.61) 


and the limit is accelerated for standard 0(a) improvement. Similar relations apply for the 
short-distance coefficients of HQET currents of dimension five and higher. 

Let us summarize the main results of this section. We have given a description of heavy- 
heavy lattice currents in HQET, which parallels the description of continuum QCD. The 
parallel structure shows that cutoff effects are isolated into the mismatch of short-distance 
coefficients for lattice gauge theory and continuum QCD, Eqs. i |2. 1 l| i and (|2.43| ). If the lattice 
currents (and Lagrangian) have enough free parameters, one can reduce the mismatch and, 
in this way, bring lattice calculations closer to continuum QCD, term by term in the heavy- 
quark expansion. In the practical region for the charmed quark, m c a ~ it is also worth 
remarking that (for Wilson fermions) one could also apply the heavy-light formalism of our 


3 One must, of course, take the same lattice currents. The improved currents of Ref. [9] and of 
Sec. Ill are not the same. 
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companion paper 0 . Then one sees, extending Eqs. (|2.56|) and (|2.57| ) to higher dimension, 
how the whole heavy-quark expansion for a quark with mass m c a < 1 is recovered. 

The remainder of this paper pursues the program of HQET matching in perturbation 
theory. One-loop corrections to the rest mass m\ and the kinetic mass m 2 have been con¬ 
sidered already in Ref. [|33|]. We construct lattice currents suitable for matching through 
dimension-five operators in the Lagrangian and dimension-four in the currents. Our heavy- 
heavy lattice currents turn out to be nearly the same as our heavy-light lattice currents. We 
then calculate the matching factors Z v and Za ± at the one-loop level, which are needed to 
fix the overall normalization of the heavy-light currents. 


III. LATTICE ACTION AND CURRENTS 


In this section our aim is to define heavy-heavy currents with Wilson fermions that are 
suited to the HQET matching formalism. Our construction is similar to that for heavy-light 
currents 0. The descriptive part of the HQET formalism applies as long as the lattice 
theory has the right particle content and obeys the heavy-quark symmetries, as it does with 
Wilson fermions. To use HQET systematically to improve lattice gauge theory, however, one 
should make choices to ensure that SCi and 5Bji [cf. Eqs. (|2.43|) and (|2.44|) 1 remain bounded 
in the infinite-mass limit. For convenience we focus on the case when v' = v. Then good 
behavior is attained by mimicking the structure of Eqs. (|2.13|) - (|2.18|) , and free parameters 
in the currents can be adjusted so that bCi and 6 Bji (approximately) vanish. We show how 
to do so in perturbation theory, obtaining B V j\ at the tree level and, in Sec. |TV]. the matching 
factors Zy {l and at the one-loop level. 

A suitable lattice Lagrangian was introduced in Ref. [£| and reviewed in the context of 
HQET matching in our companion paper on heavy-light currents []2j. We will not repeat the 
formulas here, but refer the reader to Ref. It is enough to recall that the Lagrangian has 
several couplings— r s , (, cb and eg. The coupling r s is redundant, in the technical sense, 
but the others can be tuned to match lattice gauge theory to continuum QCD; in this way 
they become non-trivial functions of the heavy quark’s bare mass m^a [£J|. In particular, one 
may define the renormalized heavy quark mass by identifying C^ at = C 2 , and one can adjust 
cb so that C x % = Cb- (If one also requires = m 2 , one can adjust both ( and the bare 
mass Pf?5|.) 

The corresponding currents for transitions from one heavy fermion to another are defined 
as follows. First define a “rotated” field fO.ffl 


= [1 + adiq'y ■ Aat] i’q, 


(3.1) 


where is the field in £0 of flavor q (q = c,b), and _Di at is the symmetric covariant difference 
operator. (The lattice Lagrangian £ 0 is given in Eq. (3.1) of Ref. [0.) The bilinears 


Vq = yjrvb, 

c *y* 75^6- 


(3.2) 

(3.3) 


have the correct quantum numbers, but not enough freedom to match at the dimension-four 
level of HQET. Thus, these currents must be improved. We focus on matching with v' = v, 
so we take 
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(3.4) 


= K 


^lat — ^0 


E b viQvi, 

*G{2,3,5,6} 

E <>*<35,, 

ie{ 1,4} 


(3.5) 


where the bji are adjustable, and the higher-dimension lattice operators are 


and 


Qv 2 = 

(3,6) 

Qv 3 = 

(3.7) 

Qv 5 = 

(3.8) 

Qv 6 = 

(3.9) 

Qai = -V^ci^pAut^b, 

(3.10) 

Qaa = +vy c p±ij^5pb- 

(3.11) 


Lattice quark fields do not satisfy Eq. (|T2|) , so i/j appears explicitly. The dimension-four 
lattice operators Q are the same as for heavy-light currents, but now only corrections 
orthogonal (parallel) to v for the vector (axial vector) current are needed. An analogous 
construction for lattice NRQCD has been given by Boyle and Davies [[32]], who found the 
same pattern of dimension-four terms. 

It is worth emphasizing the difference between Eq. Q2.37|) , the corresponding equation 
for heavy-light currents (Eq. (2.40) of Ref. |2|]), and Eq. (|3.4|). Equation (|2.37|) is a general 
HQET description of any heavy-heavy lattice current. Equation ( |3.4|) is a definition of a 
specific lattice current, namely the one used in this paper (and in calculations of B —> 
transition matrix elements). In the same vein, the Qji in Eqs. (|2.13|) - ([2.26|) are HQET 
operators, whereas the Qji in Eqs. (|3.6|) - (|3.11|) are lattice operators. Correspondingly, the 
coefficients B l jf are the output of a matching calculation: they depend on the bji, which 
must be adjusted to make 5Bji vanish. The difference between Bji and Qji (i = 1,..., 6) 
from Ref. [|] on the one hand, and Bji and Qji (i — 1,..., 14) on the other, was discussed 
above. In the former, the light(er) quark is described by a Dirac field; in the latter, the 
lighter (heavy) quark is described by an HQET field. 

The difference between the definition of the lattice currents and the description with 
HQET can be illustrated by giving the matching coefficients at the tree level. A simple 
calculation of the on-shell matrix elements {c\J^ At \b) yields 


£,lat[0] _ e _( m M +m M) a /2 

1/2 ~ e UmS y r 

*at[0] _ -(m(°J + n4° , )a/2 f [ . , [0) 5 

n ~ Umfi 

Byf ] = e -< m S+ m S)«/ 2 6{;l i = 3,6, 


(3.12) 

(3.13) 

(3.14) 

(3.15) 
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for the vector current, and 


^-ylat[0] 
U A _ L 

— e -(. m iq+ m u) a / 2 ^ 



fsla,t[0] 

n Al 

— + e -( m lq+ m ib') a / 2 j 

( 1 

+ b\ 

VSmjJJ 

pla,t[0] 

-R44 

— -. e -( m lq+ m i !)“/ 2 j 

( 1 

+ b\ 

\2mj2 


Al > 


for the axial vector current. The exponentials here contain the tree-level rest mass 

m^la = log(l + m Qh a ), 


(3.16) 

(3.17) 

(3.18) 


(3.19) 


which enters through the wave-function normalization. The easiest way to derive these 
results is to combine the Feynman rule for the current in Appendix |A| with the Feynman rules 
for external lines in Appendix C of Ref. |4[], and expand the matrix element (c(p')\ J(^ t \b(p)) 
to first order in p and p'. The zeroth order yields the C coefficients, and the first order the 
B coefficients with, for our lattice Lagrangian and currents, 


2 m 


3/j. 


C(1 + m oh a) 

m oh a(2 + m oh a) 


-d 


1 h 


(3.20) 


for each flavor h. 

One should compare Eqs. (|3.12|) -( [3T8D with Eqs. (|2.32|) — (|2.34j) . The results for (7^^ 
and Cjf 1 , together with Cyj = C^ ± = 1, show that Zyj = Z^ ± = e( m ic +m i J)°/ 2 . [Recall, in 

view of Eqs. ( |2.53| ) and ( |2.54| ), that we drop the bars from Zj] To obtain SB]jJ = 0, one 
adjusts d\ and the bji. The way to adjust d \, at the tree level, is to set equal to the 
(tree-level) heavy-quark mass. Since the rest mass plays only a trivial role in heavy-quark 
physics, the appropriate (tree-level) matching condition is , which implies 


d = ((l + mpa-Q _ r s ( 

m 0 a(2 + mod) 2(l + moa) 


(3.21) 


If di is adjusted in this way, one can then take all l}j\ = 0. 

Beyond the tree level, it is convenient to define so that spatial component of the 
degenerate-mass vector current is correctly normalized. Then Q V2 and Qy 5 would be su¬ 
perfluous for equal mass, although for unequal masses both are still needed. For the axial 
current Qai and Qaa are required even for equal masses. 

For equal-mass currents it is possible to compute Z V hh nonperturbatively for all /n^a. 
One may therefore prefer to write 


Zjcb — JZycc Zybb p Jcb 


(3.22) 


and compute only the factor p JC b in perturbation theory. This split is very useful in nu¬ 
merical calculations of matrix elements j[Rj [To|j . The strong mass dependence of the Zjs 
cancels, as do the contributions of tadpole diagrams: both enter through the self energy, 
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which is common to all Zj factors. Hence, the expansion coefficients for pj are expected to 
be small, and, as a consequence, the uncertainty in pj from truncating perturbation theory 
at fixed order is under better control than the uncertainty in the corresponding Zj. The 
one-loop calculation of Zy\ h is given in Sec. |IV B| . It has been used in Ref. || to obtain 


analogous pj factors for heavy-light currents. For both heavy-heavy and heavy-light cur¬ 
rents our calculations verify that, generically, the one-loop terms p are smaller than the 
corresponding Zj '. 


IV. ONE-LOOP RESULTS 


In this section we present results for the matching factors at the one-loop level in per¬ 
turbation theory. For the Wilson and Sheikholeslami-Wohlert (SW) actions, the one-loop 
contributions were given in Ref. JT6[] . At that stage we omitted the rotation term in the 
current, which is needed to match dimension-four currents at the tree level ||T2|,[j|. 
we complete the work started in Ref. 


T6| and report results with the rotation. 


Now 
For com¬ 
parison we also present our results without the rotation, both with and without the clover 
term. After some general remarks in Sec. |1V A| , we present numerical results for the one-loop 
terms. To illustrate the general features, we present three special cases: equal initial- and 
final-state masses in Sec. EZ3 unequal masses with a fixed mass ratio in Secs. and 

unequal masses with a fixed final-state mass in Sec. |IV D| . 

In the appendices, we present expressions for the Feynman rules and one-loop integrands, 
allowing both flavors of heavy quarks to have independent and arbitrary values of all cou¬ 
plings. A computer code for generating these results is freely available |19[ . 


A. General remarks 


We shall apply the formalism of Sec. [H] with v' — v — (z, 0 ), which applies when the 
initial and final states are at rest, or nearly at rest. The matching factors Zy^ and are 
simply ratios of the lattice and continuum radiative corrections: 


^j = 


V 1 / 2 A 

Z 2 b A J 

yl /2 

Z '2c 

cont 

7V2. 7V2 

L 2 b ly J A 2 c 

lat 


(4.1) 


where Z 2 b and Z 2c are wave-function renormalization factors of the heavy bottom and 
charmed quarks. [Recall, in view of Eqs. ( |2.53|) and (|2.54|) , that we drop the bars from Zj.\ 
The vertex function A j is the sum of one-particle irreducible three-point diagrams, in which 
one point comes from the current J (J = Vji, Aj_), and the other two from the external 
quark states. The expression relating Z 2 to the lattice self energy, for all masses and gauge 
couplings, can be found in Ref. |3^]. In view of the mass dependence of the self-energy 
function, Z 2 oc e~ mia , we write 


zj = i + y ^g 2l Z^ 


j > 


(4.2) 


1 = 1 


16 



































so that the z'f are only mildly mass dependent. (This convention follows that of Ref. ||, 
but is slightly different from the one in Ref. |33|.) By construction, the exponential mass 
dependence in pj cancels out, so we write 


OO 

PJ = 1 + J29oP [ j- ( 4 - 3 ) 

1=1 

In each of the following subsections we also present the information needed to obtain the 
BLM scale [|17],|T^]. Because the ideas behind the BLM prescription are covered in our paper 
on heavy-light currents ||, we give only the essential formulas here. Let be one of our 
one-loop terms. From the Feynman diagram, it is expressed as 

c [1] = Jd 4 kf(k) : (4.4) 

where k is the gluon’s momentum. The BLM scale q* is given through 

ln(g*a) 2 = *C [1] /C [1] . (4-5) 

where 

*(44 — j In (ka) 2 f(k). (4.6) 

Studies of perturbation theory indicate that using gy(q*) an expansion parameter usually 
leads to well-behaved series, because the BLM prescription pulls in some of the higher orders. 
Here V denotes the scheme for which the heavy-quark potential is V(q) = — CFgy(q)/q 2 - I n 
other schemes the optimal scale is slightly different, so that g|(g£) ~ 9 v(q*)- 

It is also interesting to see how q* changes under tadpole improvement. If one introduces 
the tadpole-improved matching factors 

= Zj/uo, (4.7) 


where the mean link w 0 is an y tadpole-dominated short-distance quantity, the arguments of 
Ref. |R| suggest that the perturbative series for Zj has smaller coefficients. In analogy with 
Eq. (pD we write 




>“/ 2 Z, = l + £ 9 ?Z 


(4,8) 


1=1 


where 



a = 


ln[l + m Q a/uo} 


(4.9) 


is the tadpole-improved rest mass. Then 

z'i 1 = Z'J ] - - ( -i-+-i- ) 4*1, 

2 \ 1 -f- UIq c (1 1 ?7?<q5(1/ 


(4.10) 


and because z}^ < 0 and Mq 1 ' < 0 one sees that the one-loop coefficients are reduced. 
Similarly, for computing the BLM scale 
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= * Zj^ 


+ 


1 + m 0c a 1 + m ob a 


u, 


[i] 


(4.11) 


To illustrate, we take u 0 from the average plaquette, so Uq' = —Cpj 16 and *'u[ ) 1 ' = 
-0.204049(1). 

The combinations of wave-function and vertex renormalization in Zj are gauge invariant 
and ultraviolet and infrared finite. Infrared cancellation between lattice and continuum 
integrands occurs point-by-point if the masses in continuum propagators are set equal to 
the corresponding kinetic masses. Because there are no divergences, it is straightforward to 
weight the integrands with In (ka) 2 . 

Although the main results of this section, given below, are the calculation of the full 
mass dependence of the one-loop terms, let us anticipate what to expect in limiting cases. 
When both masses vanish, we must (and do) obtain the textbook result for matching in 
the Symanzik formalism. If one mass vanishes, we must recover our results for heavy-light 
matching factors. In fact, our codes are based on the same integrand functions, so this 
“check” is automatic. 

When both masses satisfy m h a -1 , both heavy-quark propagators become static, and, 
consequently, the one-loop lattice radiative corrections vanish. Thus, the matching factors 
reduce to the continuum QCD radiative corrections. Using Eqs. Q2.29|) and (|2.30|) 


lim 


777.5a, m c a—xx) 


lim 


7775a, m c a—*oo 


^ = C F 3f(m2 c /m 2b )/16K 2 , 

'll! -* cj} 1 , = Cf [3 f(m 2 Jm 2b ) - 2] /16ir 2 . 


(4.12) 

(4.13) 


The lattice contribution to the BLM numerator also vanishes in this limit, leaving 


lim 

777 . 5 a, 777. c a—XX) 

lim *j 

777 . 5 a, 777. c a—XX) 


7 LI 

Z V|| 

41] 


Cy 1 = C F 9f(m 2c /m 2b )/32n 2 + Cy \n(m 2b am 2c a), 


[i] 


(4.14) 


_ /~i 

l a , - L/ f 


f(m 2c /m 2b )~l /167T 2 T C [ }\ \n(m 2b am 2c a), (4.15) 


where the right-hand sides are due to Ref. [36 


When one mass becomes very large, but the other is held fixed, one can still deduce the 
limiting behavior from Eqs. (|4.12| )- (|4.15| ). As the heavier mass increases, m 2b a —*■ 00 , the 
lattice Feynman diagrams become independent of m 2b a. The logarithmic part then takes 
the form 3 ln(m 2 C a). This lattice logarithm cancels only part of the corresponding logarithm 
in 3f{yn 2c /m 2b ), leaving the matching factor with 3 In (m 2b a). Then non-logarithmic con¬ 
stants depend, of course, on the fixed mass of the lighter quark. In the BLM numerators, for 
m 2c a held fixed and m 2b a —> 00 , one obtains a quadratic in ln(m 2 &a), in which the coefficient 
of ln 2 (m 2 &a) is identical to those of the continuum radiative corrections, and the coefficient 
of the single logarithm and the non-logarithmic term itself are mildly dependent on m 2c a. 

Below we plot the mass dependence as a function of m^a, because it brings out the 
asymptotic behavior for both m b a = 0 and —> 00 : when m^a C 1 , m 


1 h 


m 


2 hi 


but 

when m/,a 3> 1, ~ Thus, a plot against rn^a makes it easy to look both 

at slope and curvature in the small-mass region, and at the expected logarithms in the 
large-mass region. 


I 11 presenting results below, we show together Z^j\ *z' [ }\ 
followed by the BLM scale for Zy and Z& ± . 


J 1] 
Pj ■> 


and *p^]\ for each current, 
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B. Degenerate masses 


Figure |3| shows the full mass dependence of the matching factor Z^j for the vector current, 

with rri 2 c = m 2 b, and its BLM numerator *Zyj. Figure [| shows the same for the axial vector 
current, including the one-loop term for p A± and its BLM numerator *Pa ± ■ In this case, with 
equal masses, one has simply pa ± = Za ± /Zv„- As expected, all interpolate smoothly from 

small to large quark mass. Except for * p^ L , the curves exhibit a “knee” around m^a ~ 1-2 

(■ m 0 a ~ m 2 a ~ 2-6). For *p^\ the knee effect seems to cancel, and the function crosses over 
quickly to logarithmic behavior. 

For m c a = 774,0 = 0, our matching factors must (and do) reduce to the well-known 
massless limit, obtained long ago for csw = 1 [0 and even longer ago for csw = 0 |38|]. As 
discussed above, in the infinite mass limit (with a fixed), the heavy-quark spin symmetry of 
Wilson fermions ensures that the lattice radiative corrections vanish. The infinite mass limits 


are zB = 0. *4',' = 0, zB = 4 


— —2Cf/ 16^ 2 , and *4! — — — Cf[ 4 lnfm^a) + 


1]/167t 2 . They are shown in each panel of Figs. |3] and|4] with a dashed-dotted line. As a rule, 
the tendency to the infinite mass limit becomes obvious for m^a > 5 (mod ~ 777 . 2 a > 150). 

For smaller masses (m c a, m^a ^ 1), one expects the one-loop term p[\^ to be smaller 
than Z A , because it does not contain any tadpole diagrams. This expectation is borne out 
in comparing Fig. [l](a) and (c). For larger masses, the tadpole diagram is suppressed, cf. 
Eqs. ( |4 .1 0| ) and ([Flip , and p^\ and Z A± have the same asymptotic value. By comparing 


Za\ also. 


Fig. f|(b) and (d), one sees that *p A is smaller than zj A 

In Fig. |5| we combine the information from Figs. |3] and |J and show the BLM scale q*a for 
the vector and axial vector currents. In the region of greatest interest, ran, < 1 (■ m, 2 b < 1-5), 
we hnd q*a ~ 2.8 for the clover action. In Fig. |6] we apply tadpole improvement. In the 
region of greatest interest, mib < 1, we hnd a smaller q*a ~ 1.8-2.0 for the clover action. 
The singular behavior for Zy^ for mib ^ 1.5 arises because the denominator vanishes. If 
numerator and denominator in Eq. ( fL5p have the same (opposite) sign, then q*a —> oo 
( q*a —> 0). In such a case the BLM prescription does not make sense and should be 
mo diked 0- 

With equal masses, the matching factor Zy„ can be computed non-perturbatively, allow¬ 
ing us to test how well BLM perturbation theory works. Figure |7| plots exp(— m^a)Z Vn vs. 
m^a for several methods of calculating Z v . In perturbation theory, we truncate Eq. (|4.2|) at 
the hrst non-trivial term, and use either the bare coupling g g or the BLM prescription gy(q*)- 
We also truncate Eq. (|4.8|) at the hrst non-trivial term, which is tadpole improvement. As 
a non-perturbative check, we define 


Z* p 

VII 


(^ 1 ^ 674 ^ 1 ^) 

(H\H) 


(4.16) 


where H is a 6-havored hadron. Results for several b masses at bare gauge coupling g$ = 
6/5.9 [421 are shown in Fig. [7]. We also plot the Taylor expansion for small mass 


Z v (m 0 a) = Z v [ 1 + b v m 0 a\ , 


(4.17) 
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m 2c ~ m 2b "> 2c = m 2b 




FIG. 3. Full mass dependence of the one-loop coefficients of the matching factors of the vector 
current with equal masses: (a) Zyj and (b) *Zyj. Filled (open) symbols denote the SW (Wilson) 
action. Solid (dotted) lines connecting squares (circles) indicate the rotation is included (omitted). 



(b) 





(d) 

FIG. 4. Full mass dependence of the one-loop coefficients of the matching factors of the axial 
vector current with equal masses: (a) Z^ , (b) *Z^ , (c) , and (d) *p[^ . The symbols have 

the same meaning as in Fig. |3|: dashed-dotted lines show the infinite-mass limit. 
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(b) 

FIG. 5. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 
(b) Za ± , with m c = The symbols have the same meaning as in Figs. |j and [|. 




(b) 

FIG. 6. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 
(b) Za ± , with m c = mj, after tadpole improvement. Curves from Fig. || are shown in grey. 



FIG. 7. Comparison of full mass dependence of the matching factor Zy- . The symbols are a 
non-perturbative, hadronic calculation j42|. The dotted curve is bare perturbation theory. The 
solid curves are BLM perturbation theory; the heavy (light) line includes (omits) tadpole improve¬ 
ment. The dashed-double-dotted curve shows the first-order Taylor series, with intercept and slope 
determined non-perturbatively. 


21 
























where the notation on the right-hand side follows Ref. [0j], From the two small-mass points, 
we determine Z v = 0.7247 and by = 1.5588 (for gfi = 6/5.9 and c S w = 1-50). Figure |7] 
reveals several lessons. First, bare perturbation theory deviates from the non-perturbative 
points by several percent. Second, BLM perturbation theory is quite accurate for all masses. 
The deviations are a few percent: they are neither completely absent nor several percent. 
Third, it is less accurate to neglect the full mass dependence (using only Zy from ma —> 0) 
than to use BLM perturbation theory, already around the charmed quark mass (mi ~ 0.5). 
Finally, although Eq. ( [4.171 ) remains within a few percent of the full mass dependence for 
mi < 0.5 (cf. dashed-doublc-dotted curve), it does not have the correct large mass limit. It 
gives an explicit example of the pitfalls that arise in applying the improvement program of 
Refs. |||10] to heavy quarks. 


C. Unequal masses: fixed mass ratio 


In a typical usage of the matching factors, one wants to calculate a physical process, such 
as B —> D^lu. To monitor lattice spacing effects, one would like to repeat the calculation at 
several lattice spacings with the ratio of the initial and final heavy quark masses held fixed. 
It is therefore useful to see how the matching factors vary with m^a at fixed z = m| ( (| /mjj)]. 
For illustration we choose z = 0.256, which was used in a recent calculation 
factors for B —> Dlu. 

' [1] *TW,pW,andVw, 


13 of the form 


Figure § shows the full mass dependence of Zy], *Z\, 


held fixed at 0.256. Figure [| shows the same for , *Z 
is qualitatively similar to that with z = 1. After a knee, the expected infinite mass limit 
is again reached for mf^a > 5. The asymptotic behavior is shown with the dashed-dotted 
lines: 


4i] 


n [1] 

Pa , 


and 


with z = m 2 c /m 26 
The behavior 


U 1 ! 

Pa 


and 


P [ y\ = Z§ = +0.901CV/167T 2 , 

P [ l[ = Z [ i[ = -1.099CA/167T 2 , 

*p [ y\ = *z£ j ] = C f [ 0.124 + 1.802 ln(m 2 a)]/16vr 2 , 
= *Z l 2[ = C f [ 1.248 - 2.198 ln(m 2 a)]/167r 2 , 


(4.18) 

(4.19) 

(4.20) 

(4.21) 


using /(0.256) = 0.3003. For both currents, the crossover in *pj^ is smooth enough that no 
knee is manifest. 

At the largest masses our numerical integration (with VEGAS [^9| or BASES ]40|) dete¬ 
riorates, because ultraviolet divergences of the continuum vertex and self-energy diagrams 
do not cancel point by point (as they do when z = 1). If accurate values of the one-loop 
terms were essential, this numerical difficulty could be avoided. But it arises for values of 
rriiba where the quarks are essentially static, and one may just as well make them completely 
static. 

In Fig. [T7] we combine the information from Figs. |8] and [1] and show the BLM scale q*a 
for the vector and axial vector currents. In the region of greatest interest, < 1, we find 
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m 2c = 0.256 m 2b m 2c = 0.256 m a 




m 2c = 0.256 m 2fo 




FIG. 8. Full mass dependence of the one-loop coefficients of the matching factors of the vector 
current with z = 0.256: (a) Zyj , (b) *zj^, (c) and (d) . The symbols have the same 

meaning as in Figs. [3] and [|. 


q*a ~ 2.8 for the clover action. In Fig. 0 we apply tadpole improvement. In the region 
of greatest interest, mu, < 1, we find a smaller q*a « 1.8-2.0 for the clover action. The 
singular behavior for Z V] near mu, ~ 1.5 (and for Z& ± near mu, ~ 2.5) arises, again, because 
the denominator vanishes, so one should switch to the modified prescription of Ref. [[Il] . 

As we have mentioned above, when rrih cG 1 for both flavors, heavy quark symmetries 
emerge |l]|. One consequence is that the one-loop lattice radiative corrections vanish, as 
Figs. [3]. |4|, |8]. and ^verify numerically. There is an interesting further consequence for /y- . in 
which the unphysical e mia normalization drops out. Its approach to the infinite-mass limit 
must satisfy^ 


C 


Vii 


Pv« 


_ s~ila,t _ a 2 

— U V\\ ~ ^2 


c (2) + c (3) £ 2 


(4.22) 


4 In Sec. ||, Cy* is a generic notation for a matching coefficient in HQET. Here it denotes the 
same coefficient for the specific current JZy° c Z V bb 4 r c 7 4 T b■ 
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m. = 0.256 m_, 

2c 2b 


m. = 0.256 m n , 

2c 2b 



m 2c = 0.256 m 2b 


(b) 



m. = 0.256 m,. 

2c 2b 




FIG. 9. Full mass dependence of the one-loop coefficients of the matching factors of the axial 
vector current with z = 0.256: (a) Z^ ± , (b) *Z^ ± , (c) p^ ± , and (d) *p \The symbols have the 
same meaning as in Figs. H and [J 




FIG. 10. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 
(b) Za ± , with m 2 c = 0.2567772b. The symbols have the same meaning as in Figs. || and|9|. 
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FIG. 11. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 
(b) Za ± , with 7B2 c = 0.256m26, after tadpole improvement. Curves from Fig. H are shown in grey. 



z = 0.256, 0.3, 0.5, 0.9, 0.99. 


where the are independent of the masses, and 


^2 

S 2 


1 1 

2m2 C a 2m 2 ba ’ 

1 1 

2 m 2c a 2 m 2 bd 


(4.23) 

(4.24) 


Equation (|4.22|) is a version of Luke’s theorem H 


It follows because the mass depen¬ 
dence must be symmetric under the interchange m c mb , and because py^ must van¬ 

ish when m c = mb. Following Ref. |32|] we test whether our one-loop calculation satisfies 
the theorem for a variety of different mass ratios. Figure |l^ plots vs. S 2 for 

z = 0.256, 0.3, 0.5, 0.9, 0.99. All mass combinations lie essentially on a single curve, ver¬ 
ifying Eq. ( |1.22|) . Small deviations are possible (and noticeable in Fig. |l^), arising from 
higher orders in the small p exapnsion of the quark propagator, where coefficients other 
than 1/ (2m 2 a) appear. 
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FIG. 13. Full mass dependence of the one-loop coefficients of the matching factors of the vector 
current with m,Q C = 1: (a) , (b) *Zyj , (c) f)y J, and (d) */?[/J. The symbols have the same meaning 
as in Figs. |] and [|. 


D. Unequal masses: fixed daughter mass 


In this subsection we study the matching factors, varying only the initial mass while 
holding the final mass fixed. For illustration we set rn ac a = 1. As discussed above, one 
expects the large-mass limit to exhibit some qualitative behavior seen in the heavy-light 
matching factors in our heavy-light paper [|J. Thus, looking at this slice through the function 
of two variables can illustrate how the heavy-heavy matching factors are an extension of the 
heavy-light matching factors. 

Figure [13] shows full m^a dependence of *Z [ yj, p [ y { j, and *p [ y j], with fixed m 0c a = 1 

(m 2c a = 0.857). Figure [F| shows the same for z[y , *Z$ ± , p^\ , and *p[\\ ■ As with the 
heavy-light matching factors 0, a knee is not really visible, because soon after reaching 
m^a 2 (m 2 fea ^ 6), logarithmic behavior starts to dominate. For small m&a, one finds 
that pv n (rriba = 0,rrio c a = 1) ^ 0, but with the SW action it is very small: 0.0015017(9) 
with the rotation on the c leg and 0.0054380(10) without. 


In Fig. [15| we combine the information from Figs. [13] and [14] and show the BLM scale q*a 
for the vector and axial vector currents. In the region of greatest interest, mu < 1, we 
find q*a ~ 2.8-3.0 for the clover action (and 2.0-2.2 for the Wilson action). In Fig. |I6| we 


26 










FIG. 14. Full mass dependence of the one-loop coefficients of the matching factors of the axial 
vector current with rriQ C = 1: (a) zj-jj , (b) , (c) pj^, and (d) *p^ ± - The symbols have the 

same meaning as in Figs. H and [|. 
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FIG. 15. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 

(b) Za ± , with mo c a = 1. The symbols have the same meaning as in Figs. |l3l and [L4|. 
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FIG. 16. Full mass dependence of the BLM scale q*a for the matching factors (a) Zy^ and 
(b) Za ± , with m,Q C a = 1, after tadpole improvement. Curves from Fig. |l5| are shown in grey. 


apply tadpole improvement. In the region of greatest interest, my, < 1, we find a smaller 
q*a ~ 1.8-1.0 for the clover action (and 1.2-1.4 for the Wilson action). The singular behavior 
for Z v (Za ± ) near mu ~ 1.5 (2.2) arises, once again, because the denominator in Eq. ( |4.5|) 


vanishes, so one should switch to the modified prescription of Ref. 


E. Summary 


Let us summarize the main points of this section. The heavy-heavy matching factors are 
a function of two variables, and we have plotted several one-dimensional slices to illustrate 
their size and behavior. The dependence on the masses smoothly connects the limit where 
mha -C 1 with the infinite-mass limit(s). This is typical of radiative corrections in massive 
lattice gauge theory H, and the phenomenon should, by now, be familiar from several 


explicit examples P4|- h6| , p3 ,p!q, 

There are several checks on our calculations. Numerical results have been tested with 
two or more completely independent programs. The curves behave as expected for large 
mass. When one mass is very small in lattice units, our results smoothly connect to the 
heavy-light matching factors. When both are small we reproduce results in the literature 
for massless fermions. We are confident, therefore, that our programs compute the one-loop 
coefficients correctly for intermediate values of the masses. 

It does not seem especially useful to print tables of our results. In practice |l3|jl4|], one 
will want to analyze many different mass combinations. To facilitate this analysis, we are 
making one suite of programs freely available |19||. That will allow the user to obtain the 
coefficients to suit his or her needs. As shown in Appendix [B] we have found a way to 
incorporate the rotation into the Dirac algebra in a simple way. Since the public program 
is designed around this technique, it may be relatively easy to extend it to other lattice 
actions, such as highly improved actions. 


























V. CONCLUSIONS 


In this paper we have set up a matching procedure, based on HQET, for heavy-heavy 
currents. It is valid for all m^a and m c a , where and m c are the heavy quarks’ masses 
and a is the lattice spacing. The procedure holds to all orders in perturbation theory and, 
therefore, potentially on a non-perturbative level also. It could be applied to lattice NRQCD, 
although here it is applied to Wilson fermions. With Wilson fermions it is possible to consider 
the limit of small m^a. Then heavy-heavy HQET matching agrees with heavy-light HQET 
matching when the lighter quark satisfies m c a -C 1 , and with Symanzik matching when both 
m^a -C 1. In this way, HQET matching is a natural and attractive extension into the regime 
ma 1, which is needed for heavy-quark phenomenology on currently available computers. 

The connection to heavy-light matching has an important implication for the charmed 
quark, whose mass is not necessarily large enough for the leading 1 jm c corrections to be 
accurate enough on their own. At the same time, however, the charmed mass is small 
enough that m c a < |, on currently available lattices. In this regime the mismatches of 
short-distance coefficients of higher-dimension operators in HQET are influenced by the 
standard continuum limit (for the charmed quark), so they are easily kept down to the level 
of a s m c a and ( m c a ) 2 . Thus, more and more of the heavy-quark expansion is recovered, for 
masses where the higher-dimension terms may be significant. 

Our one-loop results for the SW action are of immediate value for lattice calculations of 
form factors for the semi-leptonic decays B — > D^lu. Indeed, our earlier one-loop results |16| 
(which omitted the “rotation” terms in the current) were used for B —> D in Ref. J13J, and 
our new results with rotation were used for B —> D* in Ref. (T^j. In particular, we have 
obtained the BLM scale q* for the renormalization factors, which reduces the uncertainty 
of one-loop calculations. Similarly, computing part of the normalization factor, namely 


ZyccZ V hh , non-perturbatively reduces the normalization uncertainty even further 

An outstanding problem at this time is the one-loop calculation of the coefficients B of 
the dimension-four terms in the HQET description. A calculation of these coefficients, and 
the subsequent adjustment of the parameters bji in the lattice currents, would reduce the 
uncertainties in (future) calculations of B — > D W matrix elements. Because heavy-quark 
symmetry forbids several power corrections (see Refs. (Ij.|T4| for details), one must, however, 
also start to consider some of the dimension-five corrections to the currents. 
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APPENDIX A: FEYNMAN RULES 


The needed propagators and vertices for quark-gluon interactions are given already in 
Ref. |33|. Here we give the additional Feynman rules induced by the rotation term of the 


heavy quark. The additional rules are easy to derive by expressing the translation operator 

= t± ll /2e ±9oaA n± IJ , / 2 (Al) 


and following the methods in Ref. ||48|| . 

There are six rules to give, with 0, 1, and 2 gluons emerging from the vertex; it is con¬ 
venient to keep rules for gluons emitted from the incoming and outgoing rotations separate. 
Let the Dirac matrix of the current be T, and let 


Then, 


R(p) = 1 + Y 7r sin (p r ), (A2) 

r 

R'(p') = 1 + id[ Y 7r sin(p'). (A3) 

r 

0-gluon = R'(p')TR(p), (A4) 

1-gluon = g 0 t a di R'(p')Tji cos (p + \k)^ (A5) 

1- gluon = g 0 t°d\ cos (p - \k) i 'y i TR(p ), (A6) 

2- gluon = gl\{t°,t h }5ijdi R'(p’)T'y i ism(p + \k + \i) l (A7) 

2-gluon" = gl\{t\ t h }5 lJ d , i i sin (p' - \k - \Z)aiTR(p) (A8) 

2-gluon 7 = glt b t a d[di cos (p 1 — |f)j 7 jT 7 8 cos(p + \k)i (A9) 


where momentum p ( p ’) is the quark momentum flowing into (out of) the vertex, and k 
and £ are gluon momentum flowing into vertex. In the two-gluon rules, there is no sum¬ 
mation over i. As in Ref. fl33 |, the matrices t a are anti-Hermitian, he., U^ = exp (^got a , 
= — CF^ikt and tr t a t b = 


+a j.a 
Z-/a_7 L ij L jk 


-S ab 
2 U ■ 


APPENDIX B: DIRAC ALGEBRA 

To compute the vertex function, there are eight diagrams to consider, depicted in Fig. [17]: 
the usual vertex diagram (with the rotation inside), Fig. |T7|(a); two diagrams with the gluon 
connected to the incoming rotation, Fig. |l7|(b) and (c); a tadpole diagram connected to the 
incoming rotation [using rule ( |A7| )1 , Fig. (T7|(d); a vertex diagram with a gluon connecting 
both rotations [using rule (|A9|) 1, Fig. [I7|(e); two diagrams with the gluon connected to the 
outgoing rotation, Fig. |17](f) and (g); and a tadpole diagram connected to the outgoing 
rotation, [using rule (|A8|) 1, Fig. [L7](h)- The two tadpole diagrams, Figs. [L7|(d) and (h), 
vanish for zero external three-momentum, because t = —k and pi = 0. 
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J 

T 

sr 

V\\ 

74 

-1 


7)75 

+3 


TABLE I. The factor sp, defined by \ Y^rlr^lr = spR 


For each non-vanishing diagram, Figs. 0(a-c, e-g), dehne the integral 


r(a-c, e-g) 

1 T 



d k 1 ^-(a-c, e-g) 

r 


(Bl) 


where k is the momentum of the gluon in the loop, and k M = 2sin(A/q i ). Let the incoming 
massive quark have couplings m 0 , r s , Cb, and Ce, and external momentum p. Similarly, 
let the outgoing massless quark have couplings m' 0 = 0, r', c ' B , and c' E , and external 
momentum p'. The internal quark lines carry momentum p + k in and p' + k out. The 
integrals / are obtained directly from the loop diagrams. Then 



1 

2 


( yM _ yM yhi _ yM ) i ( r d 

\^ J 2b cont ^2b lat ’ ^ 2c cont ^2c lat J ' / ^ cont 

d 



(B2) 


from Eq. (|4.1|) . The relation between the current J and its Dirac matrix T is contained in 
Table |lj. The expression relating f° lattice self-energy functions is in Ref. [|33[ . 


The most onerous task in evaluating the diagrams is the manipulation of the Dirac 
matrices. A convenient method is to treat each quark line separately, starting from the 
initial- or final-state spinor. Then the spinor, the propagator, and the vertices can be 
written out in 2 x 2 block diagonal form, with Pauli matrices appearing in the blocks. Once 
the Feynman rules are as complicated as in the present calculation, it is easier to manipulate 
2x2 matrices of Pauli matrices than to manipulate Dirac matrices. A special advantage of 
this organization is that the rotation bracket in Eq. (A4) merely “rotates” the rest of the 


leg. We also obtain in this way, with much less effort than in Ref. [B3 



(e) (f) (g) (h) 

FIG. 17. Feynman diagrams for calculating the vertex function. The • on each side of the <8> 
indicates the rotation. 
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A further advantage is that the vertex corrections can be expressed compactly. The 
diagram with a glnon connecting the two rotations, Fig. OM , is 


Xp C) = srd\d\ E cos 2 | k r . 

r 


(B3) 


The diagrams with a gluon going from a leg to its own rotation, Figs. [L7](b) and |17](f ), are 


'-rfb) 

-‘t 



, c 

dl D 

d'^ 

dl D' 


r 

(3-\k 2 )L' + \C'Y.KS, 


(B4) 

(B5) 


where S r = sinfc r , and the functions D , L , and K r (and primed analogs) are given in 
Appendix P The diagrams with a gluon going from a leg to the other rotation, Figs. p~7|(c) 
and [L7](g), are 


7"( c ) 

±. T 

7-(g) 

,L V 




srd[ 


c_ 

D 


(3-i k‘)n[L'] + \m')Y.Ks 2 r 

r 

(3-ife 2 )K-[i] + iK[ciy/fA 2 


(B6) 

(B7) 


where Sr is given in Table |. 

/ \ / \ 

The vertex diagram, Fig. [T7| ( a), is complicated. We hnd Zp = /DD', with numer¬ 

ator 


N^ a) = TZ[Uq]TZ[Uq] - s r n[L' 0 }n[L 0 }S 2 - C('X r . (B8) 

The part Xp comes from spatial gluon exchange: 

X r = —s r (3 - \k 2 )lZ[L'\R[L\ (B9) 

+s 2 (3 - \k)TZ[V'}TZ[V}S 2 (BIO) 

+ \{Jl[V')n[U] - s T K[L']K[<:})Y,K r S 2 r (Bll) 

r 

+i (n[u'\n[v\ - ST n\(i'\n{L\)Y,K^ (B 12 ) 

r 

+1 (n[u'\n[u] - Sr s 2 n[C}n[c}) E (Bi3) 

r 

+1(1 - 4 ) [k's 2 - 3 -£k 2 si\ n[V]n[v], (bm) 


where the last term is absent for V 4 (he., when s p = 1). The rotation enters in the “rotated” 
functions 


K[U 0 \ = U 0 + d 1 S 2 L 0 , (B15) 

E-^o] = Lq — d\Uo, (B16) 

n{U] = U + d 1 S 2 (, (B17) 

77[C] = C - d 1 U, (B18) 

1l[V] = V + d x L , (B19) 

K[L\ = L- diS 2 V, (B20) 
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and similarly for primed functions. Although the vertex diagram is not easy to write down, 
the rotation modifies it in a fairly simple way, when using the 2x2 method described above. 

We have verified that these expressions are correct by completely independent calculation 
with more common methods for the Dirac algebra. 


APPENDIX C: USEFUL FUNCTIONS 

In this appendix we list the functions appearing in Appendix [B] for the action and currents 
given in Sec. B First, let 

H = 1 + mo + \r s C, k“, (Cl) 

fi' = 1 + m' 0 + | r' s ('k . (C2) 

from now on a prime means to replace incoming coupling and momentum with corresponding 
outgoing coupling and momentum. 

When the quark propagator is rationalized it has the denominator 

D = 1 — 2/icos(/c 4 + + /i 2 + ( 2 S 2 , (C3) 

where = log(l + m 0 ). 

In this calculation, the incoming and outgoing heavy quarks both have zero three- 
momentum, so both spinors consist only of upper components. This feature is different 
from the heavy-light case | 2 ] and explains, at a low level, why the only dimension-three 
currents for heavy-heavy currents have T = 74 and 7,75 only. 

To express the useful functions compactly, it is convenient to introduce first 


U = / n- e - m i° 1+ifc4 , (C4) 

U = / n-e +m i° 1 - ifc4 , (C5) 

because these combinations appear in the other functions. Then 

U 0 = u e + m i ] ~ ik 4./ 2 _ 1(2 CE cos (l/c 4 )S' 2 , (C6) 

L 0 = ( e +m i ] ~ ik * /2 + \c E cos{\k 4 )U , (C7) 

V = C 1 + \c E sin(fc 4 ) + \c B U, (C8) 

L =—U 1 + sin(fc 4 ) + \cbC,S 2 , (C9) 

K r = r s - c B cos 2 {\k r ) = [r s - c B ) + \c B kr- ( C1 °) 


Some related functions appear when the Dirac matrix connects upper and lower components, 
as in the case of heavy-light currents [^]. They will arise for heavy-heavy matching of the 
dimension-four currents, for which matrix elements with non-zero three-momentum must be 
calculated. 
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